ON GLOBAL EXISTENCE OF CLASSICAL SOLUTIONS 
FOR THE VLASOV-POISSON SYSTEM IN CONVEX 
BOUNDED DOMAINS 

HYUNG JU HWANG, JAEWOO JUNG, AND JUAN J. L. VELAZQUEZ 

Abstract. We prove global existence of strong solutions for the Vlasov- 
Poisson system in a convex bounded domain in the plasma physics case 
assuming homogeneous Dirichlet boundary conditions for the electric 
potential and the specular reflection boundary conditions for the distri- 
bution density. 



1. Introduction 

We consider the Vlasov-Poisson system in a smooth, convex and bounded 
domain 17, with the reflection boundary condition given by 

(1) dtf + v-V^f + V^(t>-V.J = 0, for (t,x,t;) e [0,oo) X X r3, 

(2) A(j){t,x) = p{t,x) = / f{t,x,v)dv, for (t, x) G [0, oo) X 17, 

(3) / (0, X, v) = /o(x, v), for (x, v) e^xR^ 

(4) f{t, X, v) = f{t, X, V*), for X G dQ, 

where v* = v — 2(v-nx)nx, Ux is the outward unit normal vector to dQ at x G 
dQ, and we assume Q has a boundary. Here / (t, x,v) > represents the 
distribution density of electrons and <j) is the electrostatic potential, /o (x, v) 
is the prescribed initial datum, and p {t, x) is the macroscopic charge density. 
We impose the Dirichlet boundary condition for the electric potential (j), 

(5) (/>(t,x) = [fxedn. 

In this paper, we aim to understand the role of boundaries in the dynamics 
of kinetic models, in particular the Vlasov-Poisson system and to develop 
tools for boundary-value problems in kinetic models. 

For the whole space without a boundary in one and two dimensions, a 
smooth solution is known to exist globally in time |13] . |19] . For the three 
dimensional case, Batt [2], Horst [8], and Bardos-Degond [Tj proved the 
global existence of classical solutions for spherical, cylindrically symmetric. 
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and general but small initial data respectively. In the case ofn = R^ with 
arbitrary initial data, the solutions of the system ([I])- ([4]) are globally defined 
in time, as it was proved in |16j as well as in [15j using different methods. 

However, in the presence of boundaries, the mathematical theory of well- 
posedness for the solutions of the Vlasov-Poisson system becomes more com- 
plicated compared to the case of the whole space. It was proved in |6] that 
classical solutions for the Vlasov-Poisson system may not exist in general 
without the nonnegativity assumption if Q is the half-space MS]_. On the 
other hand, it was also proved in [6] that even with the nonnegativity as- 
sumption the derivatives of the solutions of cannot be uniformly 
bounded near the boundary of due to the fact that a Lipschitz estimate 
for the characteristics in terms of the initial data is not possible. 

One of the main difficulties in order to solve ([I])-®, even for short times, 
is to keep track of the evolution of the characteristic curves associated to 
([1]) which remain close during their evolution to the so-called singular set, 
defined as follows 

(6) T={{x,v) enxR^ : X e on, V eT^dn} , 

where T^dO. C is the tangent plane to 90 at the point x. 

Boundary-value problems should be treated more carefully and difficulties 
due to singularity formation at a boundary may be expected [6]. Global ex- 
istence in a half-space of solutions of (P)-® satisfying the specular reflection 
boundary condition @ and Neumann boundary condition was first proved 
by Guo (cf. [7j) by adapting a high velocity moment method in [15j. The 
proof of the global existence of solutions of ([I|)-(l3|) satisfying ^ and the 
Dirichlet boundary condition ([5]) was recently proved in ^Jj. For general 
convex bounded domains but with the Neumann boundary condition for (f), 
the global well-posedness was recently shown in |12j . 

However, a global existence theory of the Dirichlet boundary problem for 
the electric potential (j) has not been given yet and this paper is devoted 
to proving the global existence of solutions to the Vlasov-Poisson system 
([I])- ([3]) with the specular reflection boundary condition (jlj) for / and the 
Dirichlet boundary condition ([5]) for (p. 

This paper combines the methods in the papers pTj and p^j to prove 
global existence of solutions for the Vlasov-Poisson system in arbitrary 
smooth convex domains with Dirichlet boundary conditions. The analy- 
sis in [I2j allows to study problems with Neumann boundary conditions. 
This is due to the fact that an essential ingredient of the argument in [12] 
is the velocity Lemma flrst proved in [6] which shows that the characteristic 
curves associated to ([I])-© cannot approach to the so-called singular set if 
initially they are outside of it. The argument allows to generalize such type 
of velocity Lemmas to Dirichlet boundary conditions, which was obtained in 
but for the half-space case. However, simple adaptations of such tech- 
niques in the half-space problem to the case of arbitrary convex bounded 
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domains do not work. The difficulty is that, contrary to the half-space case 
[llj . we cannot have the representation formula for (p, due to the incapa- 
bility of finding an explicit form of the Green function for a general convex 
domain 0. As we will see, this problem could be settled mainly by apply- 
ing refined boundary estimates for the Laplace operator, and constructing 
relevant supersolutions. 

This paper proves global existence of solutions for the Vlasov-Poisson 
system for arbitrary smooth domains with Dirichlet boundary conditions. 
The main new contents of this paper are some technical estimates that allow 
to extend the arguments of [11] to arbitrary smooth convex domains. These 
estimates require detailed control of the newtonian potential with Dirichlet 
boundary conditions as well as some of its derivatives at points close to the 
boundary of the domain. The combination of these methods with ones in 
[12] allows to prove the stated global existence results. 

The paper is organized as follows. Preliminary notations and main result 
of the global existence will be presented in Section 2 and Section 3 is devoted 
to the Velocity lemma and the corresponding linear problem. In Section 4, 
an iterative scheme for the nonlinear problem is investigated and finally the 
global bound on a key quantity for the global existence is obtained in Section 
5. 



2. Main Result 

First we fix a point a; € and denote it as x to indicate our target point. 
Notice that our goal is to see whether a trajectory starting from the singular 
set {(x, v) G do, X R^; v -Ux = 0} propagates into the interior U x M^. So, we 
may assume that x is near the boundary dQ. Let xq be the boundary point 
closest to X. By proper rotations and translations, we may set xq = (0, 0, 0), 
X = (xi, 0, 0), and ^7 C := {{xi, X2, X3) G M^; xi > 0} so that the tangent 
plane to at xq is just {xi = 0}. 

Now, we use the local parametrization near xq to define 
X = x^i{fj,i, IJ.2) - x±n{fj,i,fi2), 

so that (/xi, ^2) is the point of dQ closest to x and n(;Ui, /X2) is the outward 
normal to dQ at xy. For this x, we represent v by 

V = V\l{fii,fJ,2) - V±n{lJ.i,fJ,2), 

where tiy (/ii, /X2) = wiui + W2U2 G T^||(^^_^2)^^ the tangential component 

of V and {ui,U2} is the basis of T^.|| (^^^^2)917 given by Ui := ^^"g^^/^^^ for 
i = l,2. 
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The system of coordinates {fJ^i, x±,wi,W2,v±) provides a more conve- 
nient representation for the set of points in the phase space x that are 
close to the singular set F defined in ([6]) as in |12j . The original equation 
([1]) takes in this new set of coordinates the following form in Lemma below, 
which is in jl2j and we skip its proof. 



Lemma 2.1. The equation (Op can be rewritten for {x,v) € + Bs 
M^, and using the set of coordinates {fJ-i, fJ.2, x±,wi,W2,v±) in the form 

(7) ^/ I ^/ I ^/ I ^ "^Z I j7 ^/ _ Q 

dt ^ 1 + kiX\ dui dx± ^ * dwi dv\_ ' 
1=1 1=1 

where 

1 + /CjX_L ^-^ 1 + kjX_i_ ^-^ 1 + kjX±^ 

j/=i j=i 

where kj are the principal curvatures, bj are the coefficients e and g from the 
second fundamental form according to the notation in [18j and F*-^ are the 
Christoffel symbols of the surface dQ. The vector E = V x<P has been written 
in the form 

(9) E = E1U1+ E2U2- Ej_n{fj,i,H2) , 
where ui,U2 are defined above. 

Remark. Notice that since the domain 0, is convex, and due to the nonneg- 
ativity assumption we have F < 0. 

To prove global existence, we need to make some necessary assumptions. 

1. Compatibility conditions for the initial data. 

In order to obtain classical solutions of (JU)-© we need to impose the fol- 
lowing compatibility conditions on the initial data /o {x, v) at the reflection 
points of dn X (cf. [6], [10]). 

(10) fo{x,v) = fo{x,v*), 

(11) [v^/o {x, V*) + V^/o {x, v)] + 2E^ (0, x) V^/o {x, v) = 0, 

where E-^ (0, x) is the decomposition of the field E (0, x) given by ([9]) and 
are the normal components to dO, of the gradients Va;, V,; respec- 
tively. 

2. Flatness condition. 

We assume that /o is constant near the singular set (cf. [7] as well as 
[10]). More precisely we will assume that /o G C^''^ satisfies the following 
flatness condition near the singular set F 

(12) /o (x, v) = constant , dist ((x, v) ,T) < 5o 
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for some 6 > small. 

We need to introduce some functional spaces for technical reasons. We 
define for fi € (0, 1) , 

llfll f \Vf{x,v)-Vf{x',v')\ \ , „^„ ^V7 Y7 ^ 

^ ^ ix,v),(x',v')€nxR^ \ \X - X \ + \V - V \ J \ I 

C]'^ ifl X = {/ e C^'^ X M^) : / compactly supported, ||/|lci.M(nxffi:^) < ' 

\V^f{t,x)-VJ {t\x') 
\x — x'\ 



cV'i''([0,T]xn) = _SUP |„_ 



I V./ (t, i;) - V./ (t', x', ^;)| + I V./ (t, x, t;) - V./ (t', x', 

x,x'ef^, t,i'e[o,T] 

+ ll/llc([0,r]xnxK3) + ll/*llc([0,T]xfixIR3) 



sup \ TiTTTI Tim + 

X — x' r + m — v'y 



We define the spaces C ([0, T] x S7) , C ([0, T] x O x M^) as the spaces of 
continuous functions bounded in the uniform norm. 



The main result: Global existence Theorem. 

The main result of this paper is the following. 

Theorem 2.2. Let fo G Cg''" {Q x M^) for some < fi < I with fo > 
and let fo satisfy f70j] - f7^] . Then there exists a unique solution f G 
^tix\) ((0' oo) X O X M3)^ g c-ijS.A (jQ^ oo) X 0) , /or some < A < /x, o/ 
i/ie Vlasov-Poisson system |iP-(Oj mt/i compact support in x and v. 

3. Velocity Lemma and Linear Problem 

Next, we introduce the evolution of characteristic curves associated to the 
Vlasov-Poisson system with the specular reflection at the boundary. 

Let E{t,x) := Vx0(i,x) be given. We define {X{s;t,x,v),V{s;t,x,v)) G 
$7 X such that for each (x, u) G $7 x E^, 

dX 

(13) 'ds^^' ^' ^ ^' 

(14) -^-(■s; a^, I') = E{s, X{s; t, x, v)) = Vx(t){s, X{s; t, x, v)), 
as 

(15) (v(t),y(t)) = (x,^). 
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as long as X € ri. The reflection boundary condition says that iiX{si;t, x, v) € 

dQ for some si G [0, T], then 

(16) 

V{sf;t,x,v) = lim^V{s;t,x,v) = (y{s^;t,x,v))* = (lim^V{s;t,x,v) 

S>Sl ^S<Sl 

Here, V* = V — 2{V ■ nx)nx where nx is the outward unit normal vector 
to dn at X. 

Before giving an explicit formulation, we consider some underlying moti- 
vations. If we rephrase the Velocity Lemma, it is equivalent to saying that 
a trajectory starting near the singular set {x± = v± = 0} remains near it in 
the future. More precisely, by using the local coordinates, we represent the 
normal component of the characteristic equations from Lemma l2. II by 

dx± dv± \ , wfbi 



v± 



EAt,x) + Y,Y 



dt dt ^ 1 + kiX± 

1=1 

where E±{t,x) is the normal component of E{t,x), /cj's are the principal 
curvatures, and 6j's are the coefficients e and g from the second fundamental 
form, according to the notations in [18]. Notice that a trajectory cannot 
escape from the singular set, provided v± < near the boundary. Roughly, 
this is true because E± < due to Hopf Lemma and ftj < by the convexity 
of ri. Finally, we define a Lyapunov function 

2 2 01, 

a{t, X, v):= — - cP{t, x) - 2^ ^^^^ x± 

1=1 * 

and confirm the stability by differentiating it along the trajectory. It may 
be possible to choose another function which is equivalent to x± + v"^, but 
the functional a makes computations simpler because of the cancellations. 

Now, we begin to show the following Lemmas which will be the crucial 
estimations required to derive our main result. Recall that Q and x are 
given in Section 2. 

Lemma 3.1. Let T > 0. Suppose that 4>{t,x) solves the following boundary 
value problem 

A(j){t,x) = p{t,x) for {t,x) £ [0,T] X n, 
(j){t, x) = if X G 

where p € C"'^([0,T] x Q) is given by 

(17) dtp + V -j = in [0, T] x Q, 

for some j € (C"'^([0,T] x il))^. Then we have 

d<i) , ~, ^ ^~ ~ ,^ 
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where C > depends only on L := diamfi and \\j\\oo- 

Proof. Let R = \x\ <^ 1. We change the variables x and y hy X = and 
Y = ^. Also, X = = (1,0,0). Let G be the Green function for the given 
domain fl. Then, by the representation formula and (|17p . it is sufficient to 
show that 

(18) \VYGiX,Y)\dY < |(1 + |logi2|). 

R 

We split the region of integration into several parts. 
Case 1. If \Y\ < 4, then we decompose G{X,Y) = G{X,Y) + W{X,Y) 
where 

1/1 1 



G{X,Y) 



4ttR\\X-Y\ \X*-Y\ 



is the Green function for the Half-space restricted to ^ x ^. Here, X* 
represents the reflection of X with respect to the plane {Xi = 0}. Clearly, 
we have 

C 



I \VYG{X,Y)\dY < 

J\Y\<4 ^ 



On the other hand, let ^ = VyVl^. Then, ^ satisfies 

R 



Ax^'(X,y) = forX,yG^, 



Firstly, if dist(y, ^) > 1, then for ah X G ^, we have 

l^(x,y)| < ~ ^, < < c, 

' ^ ' ^' - ijjx-yp - ixj'^ + i - ' 

by Taylor Theorem and the quadratic approximations of Notice that 
the constant C > can be chosen uniformly with respect to Y . Thus, 

\^{X,Y)\<C for ah X eg, 

by the maximum principle. 

Secondly, if dist(y, ^) < 1, let yo € ^ be the boundary point closest to 
y, i.e., dist(y,yo) = dist(y, ^). Then for X G ^, we get 

|^(X,y)| < ^|^^y|2 if |X-yo| >lor <Ci?, 
by the triangle inequality and the convexity of Vt. And, we have 

|^(x,y)! < — ^ ifCi?< |x-yo| <i. 
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by Taylor Theorem. Putting these together, we let 

(C 



R 



C 



R\X-Y\'^ 
C 



if |X-Fo| > 1, 
if |X-yo| <CR, 
if CR< \X -Yq\ < 1. 



Now, for |y| < 4 and dist(y, ^) < 1, by the maximum principle, we con- 
struct a supersolution via the Poisson integral formula, 



\^{X,Y)\ < 



1 



J\£.-Yo\<CR JCR<\i-Yo\<l J\i-yo\>^ 



-Yo\<CR JCR<\i-Yo\<l J\i-Yo\l 
Let r] •.= Y — Yq. Then \X — Y\ = \{X — Yq) — r]\. Using the inequality 



C 



-Yo\<CR (1 + |e|2)5 R\i^ - Yo) - ri\' 



c 



R JmcR ICP + \ri\' 



< 



c 

R 



log 



CR 



dist(y,f ) 



we obtain 



and 



MX,Y)\<^ 



log — 



CR 



dist(y, f ) 



_^ist(F,i^)<l 

|Yt<4 



anw, i^'(x,y)|dy<^(i + |iogi?|). 



Case 2. For |y| > 4, we fix a point Y = Yq. Rescale the variables by 
?7 = ^ and ,e = ^ so that ^, r/ € Define g{^,ri) := G{\Yo\^, \YQ\ri) = 

G{X, Y). Since AxG{X, Y) = S{X -Y),we have 

by a change of variables. Let <^(^,?7) := \Yo\Vr]g{^,rf)- Then we have 
A^¥'(e, rj) = V^(5(C - v) for all ^, rj e 
^{^,V) = if^G^. 
Now, we again divide into two cases. 

• If dist(ry, j^Tr) — TIJ' define ipiC^rj) satisfying 



\Yo\R^ 



1 



1 



47r ^|C-r?| 
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Since A^(/? = V,;5(^ ~ '^)) wc have 

A^i;{^,r])=0 for all ,e, ^7 € 

By assumption, \'4>{^,'n)\ ^ C for all ^ € j^^- -^^ maximum 
principle, we have 

m,v)\<C foralUGpg^. 

If we apply the boundary regularity theory to the restricted region 
G 1^1 — 5} with \ri\ = 1 fixed, then we have 

|VeV(^,r?)|<C foralUGpg^. 
Hence, we obtain 

and 

\VyG{X,Yo)\< ^ 



• If dist(ry, j^^) < we fix r/, and say rjQ. Let 

r = e + 2dist(e, {r? G M3; (r? - fJo) • i^{f]o) = 0})u{rjo), 

where ffQ is the boundary point closest to rjo, and z^(r?o) is the outward 
normal vector at r/g. This means that is the reflection of ^ with 
respect to the tangent plane at TjQ. We define w{^, 77) such that it 

satisfies 



Then we have 



A^w{^,ri)=0 foralU, r/G^, 



WoTi 



Now, for ^ G 1^^, we have 

C 1 

k(C,??o)| < 1777 iflC-%l>e' 
by the triangle inequality, and we have 

cm -Vol'' 



l'i^(^,'?o)| < 



3 
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by Taylor Theorem. Let D := dist(r7o, j^^)- We again use the 



Poisson kernel estimation to get, for ^ G B_^{0) fl j^^, 



^ f 6-^0,1 CR\z'\ 2 / 



< c 



(^ + ^|logdist (r/o,^)|). 



Here, we use the prime notation to indicate the second and third 
components in the Cartesian coordinate system. Similarly to the 
former case, we apply the boundary regularity theory and obtain 



■( 



|Vs«.(5,>»)| SCIi^^ + flllogdist ta.ifSjj) 
for all ^ e B_i^ (0) fl jy^^- Hence, by the change of variables, we get 



1^01 



|VyG'(X,yo)| <C — ^ + — ^llogdist (-^, 



In conclusion, if \Y\ > 4, we have 

/ \VYG{X,Y)\<fY 
Ja<\y\<^ 



For the second term, we get 



L 



^niogdist 



4<|y|<||yP' " '\Yy\Yo\R' 

|logfl| 

y4-2"<|y|<4-2"+i \y\ ' ' ' °' 

logfll 



KCr'y,'^'' [ Uogdist {^,^)\d^Z. 



„=0 •^4<|Z|<8 

< c 
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Notice that the last integration is bounded. Thus, we have estabUshed (fT8]l 
for \Y\ > 4. □ 



0x2 



Lemma 3.2. With the same assumptions of Lemma \3.1\ we have 

< Cxi{l + [ fog xil), 
where C > depends only on L and \\p\\oo- 

Proof. As in the proof of Lemma 13. H we take the scaled variables X = ^ 
and Y = ^ where R = \x\. Say X = ^ = (1,0,0). To compute \-§^\, we 
divide it into two cases. 

Case 1. If \Y\ > 2, we decompose the Green function G{X,Y) = 
G{X,Y) + W{X,Y) where 



G{X,Y) 



1 



1 



1 



4ttR \\X -Y\ 



\X - Y* 



and Y* represents the reflection of Y with respect to the plane {Xi = 0}. 

C 



Notice that if |X| < ||y|, then we have 



\G{X,Y)\< 



R\Y\ 



Moreover, since < W{X,Y) < -G{X,Y), we get 

G 



\G{X,Y)\< 



R\Y\ 



Now, we take the variables C = 1^ ^^id r/ = and consider the restricted 



region JIq = {'^ G -j^; \^\ < |} with Y fixed. Since 



A5G(e,r?) = Oinf)o, \G{^,v)\ < 



G 



RY 



on dQo, 



applying regularity theory leads to 

G 



< 



R\Y\ 



for any multi- index a. 



Let ,^ = Since |^(0,r/)| = 0, we have 



and 



dG 
dG 



dX-2 



{X,Y) 



< 



G 



R\Y\^' 
G 



< 



R\Y\3' 



Case 2. If |y| < 2, we denote as Yq the boundary point closest to Y. 
Decompose G{X, Y) = G{X, Y) + W{X, Y) where 

, 1/1 1 
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Here, Y is the reflection point of Y witli respect to the tangent plane at Yq. 
For X e we express W{X, Y) in terms of e± := \y-Yo\ ' ^ \y-yI\ ' ^'"■'^ 



V 



R 



as 



W{X,Y) 
Then we have 

1 



4:TtR\Y — YQ\\\e± — r]\ \e± — r]*\ 

' <CR\Y -YoWtjI'^ if|??|<l, 



1 



|e± - til 



\e± — T]^ 



< 



CR\Y-Yo\ 
1^1 



if \v\ > 1- 



Combining these, we obtain Y)\ < C. On the other hands, for X € ^ 

on the line segment \X\ = 4, we have y)| < ^. Indeed, this can be 

done by taking — G as a supersolution, and then applying the maximum 
principle. Now, if we consider the region ^ n {|X| < 4} with Y fixed, then 
using the regularity theory leads to 



\VxW{X,Y)\ < 



and by adding |VxG| term, we have 

""AY) 



dX' 



< 



C 
R' 

C 



From all these calculations, we conclude 
dG 



dX' 



iX,Y) 



dY 



+ 



\Y\>2 J\Y\<2 



i?|x-y|2 

dG 



dXo 



iX,Y) 



< 



G 



\Y\>2 



R\Y\'^ 



dY + 



dY 
G 



\Y\<2 i?|x-y|2 



dY 



G 



<-(l + |logi?|). 



□ 



We now give the main result, which plays the role of Velocity Lemma in 
our setting. 

Lemma 3.3 (Velocity Lemma). Let {X{s;t,x,v),V{s;t,x,v)) be the char- 
acteristic curves associated to the Vlasov-Poisson system defined previously. 
Suppose 4){t, x) satisfies the assumptions of Lemma I3.il Then there exist 
constants Gi and C2 > depending only on ^, \p\l°° , md ||j||L°°, such that 
if Xj_ is small enough, we have 

Ci(Xx + Vim < {X^ + Vl){s) < G2{X^ + Vim, 

fors,te[0,T]. 
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Proof. Due to Hopf Lemma, we can choose the constant eg > such that 
(19) ^{t,x) <-eoXi_ 

for x± small. We define 

a{t, x,v) = -- c^it, X) - j x^, 

where ftj's are the coefficients of the second fundamental form, and fej's 
are the principal curvatures of the surface di}. Notice that ftj < by the 
convexity of fi. So, a{t, x, v) is equivalent to x_\_ + v\. That is, it is sufficient 
to show 

Cia{t,X{t),V{t)) < a{s,X{s),V{s)) < C2ait,X{t),V{t)) 

for s,t S [0,T]. By differentiating a with respect to t along the charac- 
teristics and by representing the field E as E = Vx4> = Eiui + E2U2 — 
E±nx{lJ'i, we have 

da d<p bi 

— {t,x{t),v{t)) = - — (t,x) - 2_^WiEi — 

i * 



-Y^y^w^b^E. + w,- -2^ ^ ^ ^^^^ 

I j,i 

, Sr-f 1 d,ki\ x±wfbi 

+ Y[v^h-x^-)j^-^^, 



where Tj/SS are the Christoffel symbols. Using Lemma 13.11 13.21 and the 
equation (fT9]l . we obtain 

da 



^^it,X{t),V{t)) 



< Ca{l + |loga|). 

Therefore, the Lemma follows by Gronwall inequality. □ 



We give the theorem of well-posedness for the linear problem ([I|) , dS]) , ([!]) 
in the following theorem. 

Theorem 3.4. Assume that E G C^^ff ([0,r] x Cl) for some /i G (0, 1) . 
Suppose that fo G Cq''* [Q x M'^) for some fi > and fo>0- Then there 
exists a unique solution, f G C\.'(^^\-^ ([0)^] x x , to the linear Vlasov- 
Poisson system for some < A < /x. Moreover the function f 

satisfies 

(20) / > 

(21) f {t,x,v)dxdv = fo{x,v)dxdv , tG[0,T]. 
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Proof. The key point in the proof is that the characteristics (I13p - ()16p inter- 
sect the boundary dO, x M'^ at most a finite number of times and they never 
intersect with the singular set, which is due to Velocity lemma and Lemma 
3.31 The essential procedure is similar to the proof of Theorem 2 in jT2] , we 
skip the details of the proof. □ 



4. Iterative approach for the nonlinear problem 

In this section, we will show the global existence of classical solutions to 
the fully nonlinear Vlasov-Poisson system ([I|)-(l5]). Since the procedures are 
similar to ones in |12j . we will not try to give every detail of the proofs and 
instead we refer to [12| whenever we need. 



4.1. Iterative procedure. We will obtain a solution of the nonlinear sys- 
tem ([I])-® as the limit of a sequence of functions /" that are defined by an 
iterative procedure. More precisely, we define 

(22) f {t,x,v) = fo{x,v) , t>0,xen,veR^ 

(23) 

jn—l ( ^\ / ^n— 1, 



ft+v 




= 


(24) 




= p' 


(25) 




= 


(26) 


r (0,x,z;) 


= /c 


(27) 




= / 



(x) = / P'^dv , X € , t > 

JR3 



,x,v*) X edn , V eM'^ , t>0 

for n = 1,2, .... We assume that /o satisfies the nonnegativity condition as 
weh as (fTUD-(fT2D. 



We will use the notation 

(28) = V(/>". 

The goal is to show that the sequence converges as n ^ oo for all 
< t < CO. To this end we need to show as a first step that this sequence is 
globally defined in time for each n > 0. 

4.2. The iterative sequence {/"} is globally defined in time. Given 
a function g : Q M, we will denote as [-Jq ^^.^ the seminorm 

r 1 _ |g(x) -g{y)\ 

[9\o,X;x = sup — ■ — . 

x,yen \x — y\ 

We define 

(29) Q {t) = sup {\v\ \ {x,v) £supp f{s), < s < t} . 
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Proposition 4.1. Let /i, A G (0, 1), satisfying fi > X. Let Jq G Cq'^ x M'^) , 
/o > satisfy ^2\) . Then, the sequence of functions /" is globally de- 
fined for each x £ Q, v £ M.^ and < t < oo. Moreover we have f^ G 
Ct;'(L'i) ([°'^] X 17 X ]R3) for any T> and ||^ = ||/o|L , /Pn dx = 
f fo {t, X, v) dxdv. 

Proof. This proposition can be proved using Theorem 13.41 and by induction 
on n. We omit the details. □ 



4.3. The sequence {/"} converges to a solution of the VP system 
if the sequence {Q"} is bounded. We define the following measure for 
the maximal velocities reached for the distribution 



(30) = sup{|t7| I (x, f ) G supp /" (s) , 0<s<t}. 

Proposition 4.2. Under the assumptions of Theorem \2.Si suppose that 
Q" it) <Kforn>no, 0<t<T. Then, /" ^ f in C^^'^^^^^ ([0, T] x n x R^) 

as n— oo with < X < fi, <u <1 and where f G C'l''^^^^-^ x Q x M^) 

is a solution of 

Proof. The proof of the Proposition 14.21 is similar to that of Proposition 3 
in [12] and we omit it. □ 



4.4. Prolongability of uniform estimates for the functions 

Proposition 4.3. Let Q", Q be as in [3U\) . [2^) respectively. Suppose 
that max {sup„>„Q (t) , Q (t)} < K for < t < T. We also assume that 

f in C^.'^f^^ ([0,T] xQxR^) for any < X < < < 1. Then 

lim„_j.oo (t) = Q (t) uniformly on [0, T] . 

Proof. The proof relies on Velocity lemma, Lemma 13.31 The characteristics 
starting in a (0) > C6o remain during their evolution in the set {a (t) > CSq} 
due to Lemma 13.31 Therefore, these characteristics remain separated from 
the singular set from which we can deduce the convergence of (t) to Q (t) 
as n — >■ oo. For the details, refer to [12]. □ 



The following Proposition concerns the prolongability of the uniform es- 
timates on (t) and we skip the proof. 
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Proposition 4.4. Suppose that for some T > there exist K > and 
no > such that for any n > uq and < t < T we have (t) < K. Then, 
there exists Eq = Eq {K, ||/o||oo) > such that for < t < T + eq and n > hq 
the following estimate holds 

Q" (t) < 2K. 

We give in tlie following some of basic energy estimates for the Vlasov- 
Poisson system (cf. [4J). Proofs are standard in kinetic theory and we omit 
them. 

Proposition 4.5. Suppose that f is a solution of l^-^ defined inO < t < 
T with f{0,x,v) = fQ{x,v), where fo > 0. There exists C depending only 
on T and on the regularity norms assumed for /q in Theorem \2.2\ such that 

(31) sup / v'^f {x,t)dvdx < C 

o<t<T Jn 



(32) sup \\p{t,-)\ 

0<t<T L 



< c 



L^{n)_ 

(33) ll/(0llLp(nxiR3) = ||/o|Ilp(oxIR3) , for all I < p < oo. 

(34) ^( [ v'^fdxdv + / {Ef dx] = 
dt yJnxE.^ Jn J 



5. Global bound for Q {t) . 



In this section we show that the function Q (t) can be bounded in any 
time interval < t < T and therefore that the corresponding solutions of 
((l])-([5]) can be extended to arbitrarily long intervals. The global-in-time 
bound on Q (t) was first proved by Pfaffelmoser (cf . [I6] ) in the case of the 
whole space and the method of Pfaffelmoser has been adapted to the case 
of bounded domains with purely reflected boundary conditions at dO, (cf. 
|12j ) . The main content of the result is a uniform estimate for Q (t) as long 
as / is defined. 

From the definition (j29p of Q (t) , we obtain the following estimate 

(35) llp|loo<ll/llooQ(i)'- 
where /9 is in ([2]). 

The main result of this section is given in the following. Since the theo- 
rem and its proof do not depend on the boundary conditions for the electric 
potential (j) (whether Dirichlet or Neumann), they can be proved as in Theo- 
rem 3 in pjj . We state the theorem without its proof and without auxiliary 
lemmas. 
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Theorem 5.1. Let fo G C^'^ [Q x M^) with < fi < 1. Suppose that f G 
^i'ixl) {[0,T]xnx M3) is a solution of ^^-^ with A G (0, 1) , < T < oo. 
There exists a (T) < oo depending only onT, Q (0) , and ||/o|lci.M(r2xR3) such 
that 

(36) Q{t)<a{T) , 0<t<T . 



Proof of Theorem \2.2[ The proof is similar to that in [T2] and we omit 
it. □ 
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